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^ ■ Abstract 

The Faraday rotation angle Q is calculated in a diffusive regime in a three dimensional disordered 
slab. It is shown that tan @ is (i) an oscillating function of the magnetic field or the medium's 
internal properties, and (ii) proportional to the ratio of the inelastic mean free path /j„ to the mean 
free path I, that is to the average number of photon scatterings. The maximum rotation is achieved 
at frequencies when the photon's elastic mean free path is minimal. We have obtained the rotation 
angle of polar backscattered light taking into account the maximally crossed diagrams. The latter 
leads to an ellipticity in the backscattered wave that can serve as precursor of weak localization. 
The critical strength of magnetic field Be ^ Qc ^ A// beyond which rotation in backscattered 
wave disappears. The traversal time of an electromagnetic wave through the slab is estimated 
in a diffusive regime. The disorder enhanced the traversal time by an additional factor Im/l, in 
comparison with a free light propagation time. Comparison with the experimental data is carried 



C^ , out. 



PACS numbers: 41.20Jb,71.55Jv,71.45Gni 



I. INTRODUCTION. 



Q 



Since Anderson's pioneering work [l|, where the concept of locahzation of an electron 

moving in a random potential was introduced, propagation of electromagnetic waves through 

disordered media have been attracting a lot of interest (see, for example, Ref. 2| and 

references therein). One of the reasons for such interest is the possibility of observing 

weak localization effects in disordered systems. Most of the papers were devoted to the 

backscattering peak where weak localization effects are manifest (see, e.g., Refs. [sHTJ and 

I. 
for a recent experiment [8(). The influence of polarization effects on the backscattering 

peak was investigated in Ref. |9[. Along with the transport properties of disordered media, 

magneto-optical effects in such systems capture broader attention. Particularly, in Ref. 10|, 

the light backscattering peak in parity-non conserving disordered media was investigated. 

In such a medium rotation of the polarization plane is possible. 

A Faraday rotation is a magneto-optical phenomenon that rotates the polarization of 
light. It occurs as an internal property of medium as well as under the external magnetic 
field. In both cases the dielectric permittivity tensor of the system becomes anisotropic. 
For example, in GaAs crystal large Faraday rotation [11] origins because of the spin-orbit 
coupling. Faraday rotation is observed in astronomy. There it is used for measurement of 
magnetic fields 12|. 

In the past decade many experiment al pap ers [l3|-[l8| have been devoted to the Faraday 



rotation in disordered media. In Refs. [15Nl8| attention was drawn to the fact, that a large 
enhancement of the Faraday rotation can be obtained in nanomaterials incorporating sev- 
eral nanoparticles and their composites. A garnet composite film (mostly a Bi-substituted 
yttrium iron garnet (Bi:YIG)) incorporating gold nanoparticles or a nanoparticles solution 
with gold coated Fe20^ ISMlSl] show a change in the sign of the Faraday rotation at the 
wavelengths corresponding to the surface plasmon absorption peak of the gold nanoparticles. 
Around the local surface plasmon resonance wavelengths a resonant transmittance attenu- 
ation and a large resonant enhancement of the Faraday rotation with a narrow bandwidth 

n n 

was obtained in Ref [15|. Recently a large Faraday rotation was observed [19| in a ultrathin 
graphene film. The analogy with the Hall effect is revealed. 

The standard effective medium approximations (Maxwell-Garnet and Bruggeman the- 



ories 



20] describe the macroscopic properties of a medium using the average dielectric 



constant determined by the relative fractions of its components. These theories show only 



qualitative agreement with the experimental data of Refs. 15Nl8|. Many of the questions 
that have been raised during the experiments are still waiting for answers. For example, 
no corresponding analytical calculation has been available so far for a discrepancy in the 
bandwidth between the transmission attenuation and the resonant Faraday rotation (Ref. 



15|). Of particular importance is the influence of the disorder, induced by the presence of 
the impurities, on the Faraday angle and on several related quantities. 

In the present paper we develop a theoretical approach for investigating the Faraday 
rotation for transmitted and backscattered waves in three-dimensional disordered media. 
We show that the maximum Faraday rotation angle of the transmitted wave is achieved at 
the frequencies where the photon elastic mean free path is minimal. In the backscattered 
wave the maximally crossed diagrams contribution to the Faraday rotation angle is very 
important. It leads to an ellipticity of the backscattered wave in the localized regime. 

The paper is organized as follows. In Sec. II we formulate the problem. The coherent 
potential approximation is discussed in Sec. Ill and appropriate contribution to the Faraday 
angle is calculated. In Sec. IV we present result for the diffusional contribution to the 
Faraday rotational angle. Diffusional and maximally crossed diagrams contributions are 
taken into account for Kerr effect in Sec. V. Obtained analytical results are compared 
with experimental data in Sec. VI. Finally, the main conclusions are summarized in Sec. 
VIII. The Paper contains two Appendices, which present some technical details useful for 
the understanding of diffusional and maximally crossed diagrams contributions to the Kerr 
effect. 

II. FORMULATION OF THE PROBLEM 

Let us consider the incidence of a linearly polarized electromagnetic wave on the optically 
active medium, see Fig.l. We assume that the incident field is polarized on Oy, the plane 
of incidence is xz, wave vector is directed on Oz. Because of the scattering transmitted 
field contains all the directions of wave vector. In order to determine the Faraday angle one 
should separate the coherent part that is the one directed on Oz, see Fig.l. Faraday rotation 
angle is determined as follows 

tane=<-f>, (1) 

Ey 




FIG. 1: Geometry of the problem. Incident wave is polarized on Oy. After transmission through a 
slab the polarization vector rotates. 

where E^j-y are the components of transmitted coherent wave and < ... > means averaging 
over the reahzations of randomness. The phases of electric field components are strongly 
fluctuating in contrary to intensities. Therefore for a weakly scattering medium one can 
approximately take 

tanO 



<EyE;> 

Maxwell equation for the electric field has the form 

V^Ei{f) + £,^{f)—E,{f) = j,{f), 



VE = 



— * 

where j is the source term and the dielectric tensor is determined as follows 

eij{r) = (1 + e{f))6ij - ieijkQk- 



(2) 



(3) 



(4) 



Here Cijk is the antisymmetric tensor, g is the gyration vector directed on the magnetic field 
direction when external field is applied. In the opticaly active medium without external 
field g = fq and is parallel to the direction of photon propagation. 6{f) is the random part 
which is assumed Gaussian distributed with ^-correlation function 



< e{f)e{f') >= 75(r — r*). 



(5) 



In order to carry out averaging over the randomness it is convenient to express all the 
quantities through the Green's functions 



V2Gg(f,0- 



dr-idk 



u 



+ e^,{r) GU^ = S,,S{f-f^). 



(6) 



As was mentioned, the Faraday rotation angle in Eqs. ([T]) and ([2]) is determined by the 
coherent intensity. Hence, to separate the coherent intensities let introduce the intensity 
tensor at observation point R and in the direction of unit vector s following (9[ 

Iij{R, s) = 2 I q^dq I dfe-''^'' < Ei{R + f/2)E*{R - f/2) > . (7) 

Using Eq.(I7]) the Faraday angle, Eq.([2]), can be represented in terms of a intensity tensor 

hj{R,s) 

taue«|2^. (8) 

Here s^ is the unit vector on z and L is the thickness of the slab. Because of the symmetry 
the intensities depend only on the z coordinate. 



III. COHERENT POTENTIAL APPROXIMATION 

Averaged Green's functions take into account the fact that the right-hand and the left- 
hand circular polarized lights see different refractive indexes in the medium [lO| 






1 / (- ^ ^ 



O/cq - g^ ± iko/l 



(9) 



where fco = ^^/c, q is a unit vector on g, / = 67r/7 is the photon elastic mean free path. The 
superscripts R and A denote retarded and advanced Green's functions, respectively. 

Faraday angle consists of two main contributions. One comes from the coherent potential 
approximation and another from the diffusional scattering of photon. Coherent potential 
contribution reads 

ta„e, = ^, (10) 

where 

E,{R) = jdR'< Gt{R, R') > Jm{R'). (11) 

The source term j {r) is chosen in such a way that the coherent part of transmitted electric 
field for (? = be equal 

Ei{f^ = S^ye'^^'^'-''^'. (12) 

Using Eqs.(|9]) and f lT2|) . provided that /cq/ ^ 1, one finds 

ji{q) = 2ko{27r)H6iy6{q,)5{qy). (13) 
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FIG. 2: Dashed line is the correlation function of random field Eq.dS]) and the solid lines are the 
photon retarded and advanced averaged Green's functions Eq.(l9]). 

Substituting this expression into Eq. (fTT]l and Eq.( 1T0|) in the limit g <^ 1 and k^l ^ 1, we 
obtain 



tan 6p = e 



ikoL 



sm ■ 



ig gkoL 

— cos 

2 2 



(14) 



where g = g^. It follows from Eq. lIT^ that coherent potential approximation contribution, 
under the conditions g <ti I and k^l ^ 1, becomes inessential because contains strongly 
oscillating term exp{ikQL) when k^L ^ 1. Hence, any averaging over the wavelength of 
incident photon will make this contribution negligible. Therefore below we concentrate our 
attention on the diffusional contribution to the Faraday rotation angle. 

IV. DIFFUSIONAL CONTRIBUTION 



Using Eq.(I7]), we have 



/^(/?,s) = 2j;;^q^dqJ^,,,,..e-^''''''GtiR + f/2,p,)GUP3,R-r/2) 



X 



'IJV--5"/ -JO "i ""i JR',r',r,pi..p4 
PmnhsiPl, P2, P3, P4)G^f{P2, R' + ^ mC^iR' " ^ 1'^. PA)jf{R' + ^m3g{R' - f^, (15) 

where PmnhsiPi, P2, Ps, Pa) is the sum of ladder diagrams which was found in innumerable 
papers, see for example, |9|. We will calculate the intensity tensor in the limit (7 — t- 0. 
Therefore, it is enough to calculate the diffusion propagator for g = 0. The main contribution 
to the integrals in Eq. (TT5l) gives the term containing the diffusion pole. The corresponding 



mode has the form (see for example, Ref. 9|) 

Pmnhsi.Pl, P2, P3, Pi) = 5m'Jhs5{.Pl " P2)5(P3 - Pi)P{Pl - Ps), (16) 

where P{p) = J j0^P{K)e'^P and 

Here /„ is the inelastic mean free path of a photon in the medium and it is assumed that 

A </ <^kn,L. 

Fourier transforming Eq. flT5|) . substituting the expression P{K) into it, and using Eq.f lTT|) . 
one obtains 



dK, e'^^'G^^{Kj2 + q)Giy{K j2 - q) 



/-(i?,s;) = i27/.. <i'd<i ^ — r :;..., :;. . — - {n 



where 

The integral over K^ can be calculated exactly closing the integral contour in the complex 
Kz plane in the upper half [Z > 0). Note that the contributions from the poles Kz = i/l 



and Kz = i^^/l/Uin in the sum over polarizations cancel each other in the limit (? — !■ and 
kol ^ 1. Substituting Eq.f lT9|) into Eq. flT8|) and subsequently calculating the integrals over 
Kz and q, we obtain 

lg{R, s;) = evr^^e-^/' sin^fco^. (20) 

In deriving this expression we assume that gkol <^ 1, however kol ^ 1. The quantity 
/j^(-R, Sz) can be calculated in the same way as /^(i?, Sz). In this case the g = approx- 
imation can be used, because /^(-R, s^) is already proportional to g. Substituting Eq. (ITT]) 
into Eq.l TTS]) . one has 

lyyiR, s;) = An'e-^/K (21) 

Using the above calculated expressions for I^y{R, s^), Iyy{R, Sz) and Eq.(IH]), we finally obtain 
the desired result for the Faraday rotation angle with the diffusional contribution 

tane^^— %sin(7fcoL. (22) 

2% I 

Equation (!22|) with Eq. (^^ (see below) represent the central results of this work, tan G"^ 
presents the asymptotically exact {g — > 0, kol ^ 1) expression for the Faraday rotation in 3d 



disordered systems. Comparing the two formulae, tan0^ and tan0^, we see that the latter 
oscillates less, due to the fact that (7 — t- in the argument of the sin function. Therefore, 
the diffusion contribution to the Faraday rotation. Eg. ( 12^ . will dominate compared to 
the coherent potential approximation contribution, Eq.f lT^ . This is valid provided that 
diffusional scattering is realized in the disordered system. 

We remind that Eq.( l22|) has been derived assuming that gL/l ^ 1 however gL/X is not 
necessarily small. It is worth noticing that the ratio Im/l is the average number of scatterings 
of the photon in the medium. If the thickness L of a slab is thinner than the length y/llm 
then the diffusional trajectories are cut on the system size and the above mentioned number 



of scatterings becomes equal to L'^/P (see Refs. [2l|, |22|). In this particular case, Eq. (l22l) 
can be rewritten as 

tan e"" ^^-^ sin gkoL, (23) 

taking into account that in the diffusional regime the transmission coefficient T of a dis- 
ordered slab is of order l/L. It follows from Eqs. (l22|) and ( 123|) that the maximum of the 
Faraday angle is achieved at the frequencies where the photon elastic mean free path / is 
minimal. In the case of Eq. (12^ (L < \/Uin), it is obvious that the maximum of the Faraday 
rotation corresponds to the minimum of the transmission coefficient and in its proximity. 



This type of dependence was confirmed experimentally in Refs. 15|, |l6| . However the Fara- 
day rotation peak will be steeper because of stronger dependence on l{u) (see below for more 
details). Another pecularity of our result, Eq.( !22|) . is that it implies an oscillating depen- 
dence of the rotation on g, that is on a magnetic field, or the medium's internal properties. 
Note that the non-monotonous dependence of the Faraday rotation on the magnetic field is 
observed in Ref. |l7|. 

V. ROTATION IN THE BACKSCATTERED LIGHT: KERR EFFECT 

For completeness, we now turn to the discussion of the Faraday rotation angle for the 
backscattered wave using the procedure outlined above. The main difference now is that 



one should take into account the contribution of the maximally-crossed diagrams (Ref. |23j ) 
in Fig.2, which goes beyond the ladder diagram approximation, (i.e., regular diffusion). 
Another important point to mention also is that in the Kerr effect the boundary conditions 
at 2 = are important. 

8 



For a brief sketch of the derivation of the backscattering peak note, that the diffusion 
propagator P should fulfill the boundary condition at 2; = 0. Therefore, instead of the 
former Pmnhs, "we, based on the method of images, use a modified P(p, 21,^2) = P{p,zi — 
Z2) — P{p, Zi + Z2)- Here p is a two dimensional vector in the xy plane and we assume that the 
following inequalities are met: (7 <^ 1, /cq/ ^ 1 and gk^l ^ 1. The details of the calculations 
of the backscattering rotation is given in Appendices A and B. The final result {Un = C)0, 
i.e. no absorption) reads 

tanef = --gkol tan 6^ = ^, (24) 

where Bf, ' are the diffusional and maximally crossed diagrams contribution to the Faraday 
rotation angle of the backscattered light, respectively. Note that maximally-crossed diagrams 
leads to a pure imaginary Faraday angle which indicates the onset of ellipticity for the 
backscattered wave. As is seen from Eg. (1241) the diffusional contribution is dominant in 
the diffusion regime kol ^ 1 in contrary to the maximally crossed diagrams contribution 
which becomes relatively large in the localized regime kol ~ 1. In this regime, kol ~ 1, 
the diffusional and maximally crossed diagrams contributions are of the same order and 
maximum of rotation angle is achieved. Physically, such a behavior can be understood 



within the existing re^ 



light in a medium 24 



ation between the Faraday rotation angle and the time spent by the 



25| (see below). Indeed, in the diffusive regime the backscattering 
time is larger for low scatterer concentration (or for larger photon mean free path). With 
increasing the scatterer concentration (the photon mean free path becomes smaller) the 
backscattering time increases because of the localization effects (or because of the ellipticity 
of the polarization of the reflected wave). This transparent picture serves as a basis for 
qualitative understanding of asymptotic behavior of backscattered wave in two regimes. 
Thus, the onset of ellipticity in the backscattered wave can serve as a precursor of light 
weak localization. 

Note, that in the limit gkol ^ 1 the diffusional and maximally crossed diagrams con- 
tributions to the rotation angle of backscattering wave are very small (see Eqs. (]A12p and 



TgI)). Hence for the backscattered wave the maximal rotation is achieved for g ~ 1/kol. 



In concluding this section, let us estimate the traversal time of an electromagnetic wave 



in 3d disordered media. According to Ref. 2^ the interaction time in a slab is closely 



related to the Faraday rotation. Hence, one might expect, that the traversal time, r, of an 



electromagnetic wave through a 3d disordered slab will be larger compared with the time, 



To = Lrio/c, spent in a dielectric with open boundaries. Following Ref. 2J] and expanding 
Eq. (122 p in powers of g, and keeping linear terms only we get an estimated value for time 
r ^ {lin/^)To- Thus, the disorder enhanced the traversal time r by an additional factor /«„//, 
in comparison with a free light propagation time Tq. 

VI. COMPARISON WITH THE EXPERIMENT 



Most of the experiments on Faraday rotation in disordered media 13|- 17|] are carried 
out with the films of thicknesses much less than the light wavelengths. Our approach can- 
not be directly applied to these systems because the randomness in this case is nearly 
two-dimensional contrary to our 3(i randomness. We will compare our results with the ex- 
perimental data of the work Ref. [18|], where the sample has "id randomness. First let us 
estimate the dimensionless constant g for a water noting that in Ref. 18|] nanoparticles are 
randomly embedded in the water. It follows from ( TT^ that g = {e^ — eJ)/2 = n{n^ — n_) 
and n = {n^ + n^)/2. e± and n± are dielectric constants and refraction indexes of right-hand 
and left-hand circularly polarized light at 632nm wavelength, respectively. The difference 
An = n^ — n_ is related to the Verdet constant u: 6 = uBL = irAnL/X. Taking B = IT, 
one has An = Xu/n. For a water at A = Q32nm,n = 1.33 and u = 3Arad/Tm, see for 



example, 26|. Therefore g = nAn ~ 10~^. We see that for a water at -B = IT, g is ex- 
tremely small and the limit g -^ in our theory is justified. The elastic photon mean free 
path / in a system of randomly embedded in a dielectric host medium of metallic spheres 
can be estimated as / = 87rc^/(9f/a;^), where v is the volume of a single sphere and / is 
the volume fraction of spheres in the medium |27|. We take v ~ Ana^/S, where a ^ 85nm 

n 

is the average size of particles coated with gold in an aggregate state and / = 0.035 |18| . 
Substituting these values and A = 632n,??i into expression of /, one gets I ^ 5X. Now let us 
estimate the inelastic mean free path /j„, which is defined as kn = X/{27rImeeff) |27|, |28| 
{ImSeff = 9fe2/{£i — 1)^ is the imaginary part of the Maxwell- Garnett effective dielectric 
constant and ei, 62 are the real and imaginary parts of the dielectric constant of gold at 



incident wavelength A = 632nm, respectively). Taking ei ~ —10.7 and 62 ~ 3.23, Ref. |29|. 
one has /j„ ~ 21. 4A and it is easy to convince oneself that the conditions for photon diffusion 
A ^ /(A) <^ kni^) ^ L are reahzed in the experiment [l8|. In this case the additional large 
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multiplier Im/l ~ 4 arises in the Verdet constant calculated within the effective Maxwell- 
Garnett theory. This additional factor could explain the discrepancy between the theoretical 



calculations and the experimental data, mentioned in Ref. [18|. Multiplying the Verdet con- 
stant of water 3.4 by /j„// ~ 4 and subtracting the water value, for the Verdet constant 
of solution, one gets v ~ 10.2 which is in a very good agreement with the corresponding 



experimental value |18 |. 



VII. SUMMARY. 

We have investigated the Faraday rotation in a disordered medium. The rotation angle 
is calculated to leading order of parameters gkoL, 1/kol. It is shown that in the diffusive 
regime Faraday rotation is proportional to the ratio of the inelastic mean free path l^ to the 
mean free path /, that is to the average number of photon scatterings in the medium. For 
the Faraday angle the magnetic field oscillating dependence is predicted. If the thickness L 
of a slab is thinner than the length y/llm then the Faraday rotation is inverse proportional to 
the square of the transmission coefficient. The maximum rotation is achieved at frequencies 
where the photon elastic mean free path is minimal. We have calculated the Faraday rotation 
of backscattered waves. Maximally-crossed diagrams lead to an ellipticity of backscattered 
wave. In the localized regime k^l ~ 1 the diffusional and maximally crossed diagrams 
contributions Eq. i^^ to the rotation angle are of the same order. Therefore appearing of 
ellipticity in the backscattered wave can serve as a precursor of a weak localization of the 
light. In the backscattered wave maximal rotation is achieved for g ~ l/Zcg/. Beyond this 
critical value rotation angle rapidly disappears. The traversal time of an electromagnetic 
wave through the slab is estimated in a diffusive regime. Comparison with the experimental 
data is carried out. To summarize, although the main result for Faraday rotation, Eq. 
is approximate, it nevertheless predicts correctly many of the peculiar features of most 



experimental results [l3|-[l8|, discussed in this paper. 
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Appendix A:. DIFFUSIONAL CONTRIBUTION TO THE KERR EFFECT 

We derive first the diffusional contribution to the Kerr effect, i.e., tanG^ = —'^gk^l of 
Eq. fl2^ . To do this, let us note that the diffusional contribution to the intensity tensor, 
see Fig. 2, takes the form 

/•OO f I* 

igiR^s) = 2 J q'dqj dre-'^'' j rfndf2G'f^(^+r72,fOG;^,(n, ^-r72)P(n,r2)^,(f2)^*(f2), 

(Al) 
where 

P(ri, ra) = P{p, z^ - z^) - Pip, z^ + z^). (A2) 

Here p is the projection of ri — r2 in the xy plane. P{f) is determined by three dimensional 
Fourier transform of Eq. (IT711 and the coherent field Eh{f2) is determined by Eq. (IT2i) . Going 
through Fourier transforms in Eq.f lAip . one gets 

/5'(P,s) = 2j^q'dqJ^^^5{q,^ - Kj2 - qs.)Gg,{qu)Gi^iK^ - gi,)e^^^^/(i^,), 

(A3) 
where 



'iKzZi—Z2/l 

Jo 
and 



I{K,)= I dzidz2PiO, zi, Z2)e~''^^''-'^/^ (A4) 





P(0, zi, Z2) = I dpPip, zi, Z2). (A5) 

Here we assume that Sp = 0. Substituting Eqs.(1T71) and ( 1A2I) into Eq.l JASp . we obtain 

P(0, Zi, Z2) = ^[ZI + Z2- \Zi -Z2\]. (A6) 

Note that the expression ( ]A6p was derived ignoring the absorption and assuming that Im = 
OO. Substituting Eq. (lA6l) into Eq.( ]A4l) . one has 

„;.,) . !^,.-«« - 1) . ^. (A7) 

Note, that in order to arrive to Eq. (1A7I) . we cut the diverging integral on the upper limit at 
the system thickness L. For the backscattered direction s^ = — 1 and for R = the intensity 
I^y{R, s), Eq.f lASp . can be written in the form 

lg{0,-l)=2J^ q'dqj ±-±GUKj2-q)Giy{Kj2 + q)I{K;). (A8) 

12 



Substituting the Green's functions, Eq.(|9]), into Eq. OASP and keeping only linear terms on 
Kz, one has 



where 



jfm,-i) = /((?)- /(-^), (A9) 



oo 



We will integrate Eq. (lAlOp over K^ closing integral contour in the bottom half of complex 
plane where the function I{Kz) has no poles. The contributions to the integral come from 



the poles K^ = -iiz^Mz^! _ iM and K^ = -^IziizslM _ IM. Using Eqs.dM]) and (jM]) and 



q ql ^ q ql 

taking integral over K^ in the limit ^f — > 0, we find 

j^ lv\ — 1 ^ = —IMlM r°° 2!^2 |_ 

24£7fcg TOO 2!^g ('Alii 

I Jo [(g2-fc2(l-g))2+4fc2//2j [(g2-fc2(l+g))2+4fc2/i2j " ^ ^ 

Taking in Eq.l JAlll) the remaining integrals over q in the limit kol 3> 1 we finally obtain 

tD,. _^^ ^ 97r^lgkol 

In the final step using Eq.® and dividing /^(O, —1) by backscattered intensity at (7 = 0, 



B 



-^ro(0; ~1) = 3717/ [9|], we arrive at tanG^ of Eq. (12^ in the limit gf/co/ <^ 1. 

Appendix B:. MAXIMALLY CROSSED DIAGRAMS CONTRIBUTION TO THE 
KERR EFFECT 

In this section we derive the maximally crossed diagrams contribution to the Kerr effect. 
i.e., tan 9^ = ^ of Eq. ( 121|) . To this end, let us represent the maximally crossed diagrams 
contribution, see Fig.2, to the intensity tensor in the form 

/g(^, s) = 2 j^ q'dq U^^^.^^^...?^ e-^^^GtiR + r/2, p^)G^^ip„ R - r/2) 

PmnhsiPU P2, PS, Pi)G^l{P3, n)G^,(f2, Pi) j I (fl) j } {f2) ■ (Bl) 

It follows from the symmetry, see Fig.2, that 

PmnhsiPl^ P2, P3, Pa) = Pmshnipl, p4, p3, P2)- (B2) 
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Note that throughout the paper we consider the maximally crossed and diffusion propagators 
P^ , P at g = 0. The reason is that in the weak scattering limit k^l ^ 1 the main contribution 
gives the pole term which is unaffected by Faraday rotation or optical activity [10,]. Using 
Eqs.dlS]) and ^^, one has from Eq.dHI]) 

lf^{R, s} = 2l q^dq I dre-^'i'n dndf2Gf,iR+f/2, n)G^,ir2, R-f/2)P{n, f2)Mr2)E:{n) 

(B3) 
Going over to the two and three dimensional Fourier transforms in Eq. (IB3l) and using 
Eq.f lT^ . we have 

4^.(0,^1 =^j^q'dqj^^ '^G^^{qs„p,)G^^i-qs„p^ - 2qs^) 
J^ dzidz2Pi\qSp\, Zi, Z2) exp [-izi{p^ + ko) + iz2{kQ - p^ + 2qs^) - ^^] . (B4) 

We remind that Sp is the projection of the unit vector s on the xy plane. In the backscattered 
direction Sp = and Sz = —I- Integrating the equation ( JB4I) consequentially over zi and Z2, 

we obtain 

"+00 



where 



q'dqj -^Gg{pz)G^^{pz + 2q)F{pz), (B5) 

^^^'^ " PiPz + ko- i/2l){ko -Pz + 2qsz + i/2l){-2pz + 2gs, + i/l) " ^^^^ 



Substituting Eqs. (|9]) and f lB6|) into Eq.( ]B5|) . one has 



where 

37 r°° rfi^n r dp_ 



Jf ro, -1) = h{g) - h{-g) + h{g) - hig), (B7) 



h{9) = -1^ ir q'dq f 



2lT 



{Vz+ko-i/2l){p^-ko+2q-i/2l){p:,+2q-i/2l)(p^-yJ{l+g)kl-^iko/l) 



y 1 f RS'l 

{p. + y/{l+g)kl+iko/l){p,+2q-y/(l-g)kl-ikajl){p,+2q+y/{l-g)kl-iko/l)'' 



and 



4(^)=-#r?^^^/t 



{p,+ko-i/2l){p,-ko+2q-i/2l){p,+2q~i/2l)(p^-y/{l+g)k^+iko/l) 



X 



1 



(B9) 

{p, + y/{l+9)k^+iko/l)ip,+2q-y/{l+g)k^-iko/l)ip,+2q+y/{l+g)k^,-iko/l) 

The square roots from the complex numbers are implied in the sense of arithmetic root. In 
the next step we expand these roots on g and 1/kQl up to the linear terms 



h{9) = -H^ lo^ q'dq f '^ 



ipz+ko-i/2l)(p^-kc+2q-~i/2l)(p^+2q~i/2l)(p^-ko-gko/2-i/2l) 



X 



1 



{pz+ko+gko/2+i/2l){pz+2q-ko+gko/2+i/2l){p:,+2q+ko-gko/2-i/2l)- 
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(BIO) 



and 



^9) = -§ ir <l'dq f '^ 



X 



{pz+ko-i/2l){p^~ko+2q-i/2l){p^+2q-i/2l){p;,-ko-gko/2-i/2l) 
^ ip,+kQ+gko/2+i/2l)(pz+2q-kQ-gko/2+i/2l)ipz+2q+ko+gko/2-i/2l)- \^^^) 

We calculate the integral over p^ closing the integral contour in the bottom half of a complex 
plane. The contribution for Ii{g) comes from the residues of the following two poles 

(i).p, = -ko-gko/2-i/2l {ii).p, = -2q + ko + gko/2 - i/2l. (B12) 

The result reads 

J / N 3^7 r°° q + kp 

'^^' 8klP{gko + 2i/l)Jo \q-ko + gko/4 + z/2l){q-ko-gko/4-z/2iy ^ ' 

Eq.f lB13l) was derived taking into account the fact, that the main contribution to the integral 
on g gives the values that are close to fco- Calculating the integral Eq. (IB13P in the pole 
approximation which is justified provided that A;o/ ^ 1, one has 

''^'^ = -WH9h+m^- ^^''^ 

Correspondingly 

12iTT'~flg 

(4 + g'k^P 

In analogous manner one can show that l2{g) — hi—g) = 0. Therefore using Eqs.( ]A9l) and 
(IBlSp . one finds 



h(s)-h{-9) = ,/::':" (bis) 



Similarly to what was done in Appendix A, where the tanGf was calculated, one readily 
sees that dividing /^(O, —1) by backscattered intensity &t g = 0, lyy {0,-1) = Stt'jI J9[ yields 
tanBf of Eq. (^ in the limit gkol < 1. 
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